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m : ABSTRACT 

■ We calculate the maximum fraction of matter which is able to condense out of 

> ■ 

ly^ ' the expanding background universe by gravitational instability - the asymp- 



totic collapsed fraction - for any universe which is unbound and, hence, will 



\ expand forever. We solve this problem by application of a simple, pressure-free, 

On ■ spherically symmetric, nonlinear model for the growth of density fluctuations 

o\ 

' in the universe. This model includes general kinds of Friedmann universes, 

Oh' such as the open, matter-dominated universe and those in which there is an 

O \ extra, uniform background component of energy-density (e.g. the cosmologi- 

c/3 . cal constant or so-called "quintessence" ), perturbed by Gaussian random noise 

matter-density fluctuations. These background universes all have the property 
' that matter-domination eventually gives way either to curvature-domination 

^ \ or domination by the positive energy density of the additional background 

component. When this happens, gravitational instability is suppressed and, 
with it, so is the growth of the collapsed fraction. 

Our results serve to identify a limitation of the well-known Press-Schechter 
approximation for the time-dependent mass function of cosmological struc- 
ture formation. In the latter approximation, the mass function determined 
from the predicted collapse of positive density fluctuations is multiplied by an 
ad hoc correction factor of 2 based upon an assumption that every positive 
density fluctuation which is fated to collapse will simultaneously accrete an 
equal share of additional matter from nearby regions of compensating nega- 
tive density fluctuation. The model presented here explicitly determines the 
actual value of the factor by which any positive density fluctuation which ever 
collapses will asymptotically increase its mass by accreting from a compensat- 
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ing undcrdcnsity which surrounds it. We show that, while the famous factor 
of 2 adopted by the Press-Schechter approximation is correct for an Einstein- 
de Sitter universe, it is not correct when the "freeze-out" of fluctuation growth 
inherent in the more general class of background universes described above 
occurs. When "freeze-out" occurs, the correction factor reduces to unity and 
the standard Press-Schechter approximation must overestimate the collapsed 
fraction. 

To illustrate this effect, we apply our model to currently viable versions 
of the Cold Dark Matter (CDM) model for structure formation, with pri- 
mordial density fluctuations in accordance with data on cosmic microwave 
background anisotropy from the COBE satellite DMR experiment. For Hq = 
70kms~^Mpc^^ and matter-density parameter Oq = 0.3, the open, matter- 
dominated CDM model and the flat CDM model with nonzero cosmological 
constant yield asymptotic collapsed fractions on the galaxy cluster mass-scale 
of IO^^Mq and above of 0.0361 and 0.0562, respectively, only 55% of the values 
determined by the Press-Schechter approximation. These results have impli- 
cations for the use of the latter approximation to compare the observed space 
density of X-ray clusters today with that predicted by cosmological models. 

Key words: cosmology: theory — galaxies: clusters: general — galaxies: 
formation — gravitation — large-scale structure of the universe 

1 INTRODUCTION 

Explaining the origin and evolution of galaxies and large-scale structure and determining 
the fundamental properties of the background universe are the primary goals of modern 
cosmology. The most common assumption is that the structure we observe today (density 
structures such as galaxies, clusters, and voids, as well as velocity structures such as the 
Virgocentric Infall or that associated with the Great Attractor), results from the growth, 
by gravitational instability, of small-amplitude, primordial density fluctuations present in 
the universe at early times. These fluctuations are normally assumed to originate from a 
Gaussian random process. In this case, they can be described as a superposition of plane- wave 
density fluctuations with random phases. One important property of these initial conditions 
is that overdense and underdense regions occupy equal volumes (in other words, their flUing 
factors are 1/2). Since the density is nearly uniform at early times, overdense and underdense 
regions also contain the same mass. 

The gravitational instability scenario makes the following predictions: overdense regions, 
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The asymptotic collapsed fraction in an eternal universe 3 
because of their larger gravitational field, will decelerate faster than the background uni- 
verse, resulting in an increase of their density contrast relative to the background. If this 
deceleration is large enough, these regions will turn back and recoUapse on themselves, re- 
sulting in the formation of positive density structures such as galaxies and clusters. The 
opposite phenomenon occurs in underdense regions. These regions decelerate more slowly 
than the background universe, thus getting more underdense, and eventually become the 
cosmic voids we observe today. 

In this paper, we investigate the asymptotic collapsed fraction, defined as the fraction of 
the matter in the universe that will eventually end up inside collapsed objects. Obviously, 
this makes sense only in an unbound universe. Naively, we might think that the asymptotic 
collapsed fraction will be equal to 1/2, since half the matter is located in overdense regions at 
early times. This ignores two important effects. First, some overdense regions might be un- 
bound, and second, matter located inside underdense regions could be accreted by collapsed 
objects. The importance of these effects depends upon the particular background universe 
in which these structures form. Consider, for instance, an Einstein-de Sitter universe. In 
this case, the background density is exactly equal to the critical density, and therefore all 
overdense regions are bound, and will eventually collapse. Furthermore, it can easily be 
shown that any mass element located inside an underdense region is gravitationally bound 
to at least one overdense region. Consequently, all the matter inside underdense regions will 
eventually be accreted by collapsed objects, and the asymptotic collapsed fraction is unity. 
This is not true, however, for a background universe with mean density below that of an 
Einstein-de Sitter universe. 

Interest in models of the background universe in which the matter density is less than 
the critical value for a flat, matter-dominated universe is now particularly strong, on the 
basis of several lines of evidence which can be reconciled most economically if Qq < 1, 
where Qq is the present mean matter density in units of the critical value. (For reviews 
and references, see, e.g., Ostriker & Steinhardt 1995; Turner 1998; Krauss 1998; Bahcall 
1999). Arguments in favor of a flat universe with Qq < 1 in which a nonzero cosmological 
constant makes up the difference between the matter density and the critical density have 
been signiflcantly strengthened recently by measurements of the redshifts and distances of 
Type la SNe, which are best explained if the universe is expanding at an accelerating rate, 
consistent with Qq = 0.3 and Aq = 0.7, where Aq is the vacuum energy density in units of the 
critical density at present (Garnavich et al. 1998a; Perlmutter et al. 1998). When combined 
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with measurements of the angular power spectrum of the cosmic microwave background 
(CMB) anisotropy, these Type la SN results can be used to restrict further the range of 
models for the mass-energy content of the universe. In particular, while the SN data alone 
are better fit by a flat model with Qq < 1 and a positive cosmological constant than by an 
open, matter-dominated model with no cosmological constant (e.g. Perlmutter et al. 1998), 
the combined information from Type la SNe and the CMB significantly strengthens the 
case for a flat model with cosmological constant over that for an open, matter-dominated 
model (e.g. Garnavich et al. 1998b). Exotic alternatives to the well-known cosmological 
constant which might also contribute positively to the total cosmic energy density and 
thereby similarly affect the mean expansion rate have also been discussed, sometimes referred 
to as "quintessence" models (e.g. Turner & White 1997; Caldwell, Dave, & Steinhardt 1998). 
Such models can also explain the presently accelerating expansion rate indicated by the Type 
la SNe, while satisfying several other constraints which suggest that < 1- The results 
from Type I SNc and CMB anisotropy combined can be used to constrain the range of 
equations of state allowed for this other component of energy density px, with pressure 
Px — '^xPxC^- The current results favor a fiat universe with VIq <1 and an equation of state 
for the second component with a value oi ~ —1 (where = — 1 for a cosmological 
constant) favored over larger values of (such as would describe topological defects like 
domain walls, strings, or textures), although the restriction of the range allowed for Wx is 
not yet very precise (Garnavich et al. 1998b). 

Consider now an unbound universe with a matter density parameter Q with present 
value Qq < 1- In such a universe, the critical density exceeds the mean density, and therefore 
some overdense regions are unbound. The asymptotic collapsed fraction could still be unity 
if all the matter in overdense, unbound regions plus all the matter in underdense regions is 
accreted. This will never be the case, however. In such a universe, the density parameter VL 
is near unity at early times, and structures can grow. Eventually VL drops significantly below 
unity, and a phenomenon known as "freeze-out" occurs. In this regime, density fluctuations 
do not grow unless their density is already significantly larger than the background density. 
After freeze-out, accretion by collapsed objects will be very slow, and most of the unaccreted 
matter will remain unaccreted. The asymptotic collapsed fraction will therefore be less than 
unity. 

The asymptotic collapsed fraction is a quantity which is relevant to modern attempts 
to interpret observations of cosmic structure in at least two ways. For one, anthropic rea- 
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The asymptotic collapsed fraction in an eternal universe 5 
soiling can be used to calculate a probability distribution for the observed values of some 
fundamental property of the universe, such as the cosmological constant, in models in which 
that property takes a variety of values with varying probabilities (Efstathiou 1995; Vilenkin 
1995; Weinberg 1996; Martel, Shapiro, & Weinberg 1998, hereafter MSW). Examples of 
such models include those in which a state vector is derived for the universe which is a su- 
perposition of terms with different values of the fundamental property (e.g. Hawking 1983, 
1984; Coleman 1988) and chaotic inflation in which the observed big bang is just one of 
an inflnite number of expanding regions in each of which the fundamental property takes a 
different value (Linde 1986, 1987, 1988). In models like these, the probability of observing 
any particular value of the property is conditioned by the existence of observers in those 
"subuniverses" in which the property takes that value. This probability is proportional to 
the fraction of matter which is destined to condense out of the background into mass concen- 
trations large enough to form observers - i.e. the asymptotic collapsed fraction for collapse 
into objects of this mass or greater. MSW used this approach to offer a possible resolution 
of the infamous "cosmological constant problem," one of the most serious crises of quantum 
cosmology. Estimates of the size of a relic vacuum energy density pv from quantum fluctu- 
ations in the early universe suggest a value which is many orders of magnitude larger than 
the cosmic mass density today, and no cancellation mechanism has yet been identifled which 
would reduce this to zero, let alone one so flnely tuned as to leave the small but nonzero 
value suggested by recent astronomical observations (i.e. where the net pv is the sum of a 
contribution from quantum fluctuations and a term A/SttG, where A is the cosmological con- 
stant which appears in Einstein's held equations) (Weinberg 1989; Carroll, Press, & Turner 
1992). MSW calculated the relative likehhood of observing any given value of pv within 
the context of the flat CDM model with nonzero cosmological constant, with the amplitude 
and shape of the primordial power spectrum in accordance with current data on the CMB 
anisotropy. Underlying this calculation was the notion that values of pv which are large are 
unlikely to be observed since such values of pv tend to suppress gravitational instability and 
prevent galaxy formation. MSW found that a small, positive cosmological constant in the 
range suggested by astronomical evidence is actually a reasonably likely value to observe, 
even if the a priori probability distribution that a given subuniverse has some value of the 
cosmological constant does not favor such small values. Similar reasoning can, in principle, 
be used to assess the probability of our observing some range of values for other properties 
of the universe, too, in the absence of a theory which uniquely determines their values (e.g. 
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the value of Qq; Garriga, Tanaka, & Vilenkin 1998). In such calculations, the asymptotic 
collapsed fraction is a fundamental ingredient. 

Aside from its importance in anthropic probability calculations like these, in which one 
needs to know the state of the universe in the infinite future, the asymptotic collapsed 
fraction is also relevant as an approximation to the present universe, for the following reason. 
In an Einstein-de Sitter universe, in which there is no freeze-out, the asymptotic collapsed 
fraction is unity. In any other unbound universe, there will be a freeze-out at some epoch. If 
we live in such a universe, the freeze-out epoch could be either in the future or in the past. 
However, if recent attempts to reconcile a number of the observed properties of our universe 
with theoretical models of the background universe and of structure formation by invoking 
an unbound universe with Qq < 1 are correct, then the freeze-out epoch is much more likely 
to be in the past. If it were in the future, then the matter density parameter today would still 
be close to unity, e.g. Qq > 0.9 or Qq > 0.99.Qlf so, then the observable consequences of the 
eventual departure of the background model from Einstein-de Sitter would be largely in the 
future, as well. As such, the strong motivation for considering models with fio < 1 in order 
to explain a number of the observed properties of our universe as described above would 
vanish. In short, the current interest in a universe with f2o < 1 is consistent with a value of Qq 
small enough that the epoch of freeze-out is largely in the past. In that case, the asymptotic 
collapsed fraction should be a good approximation to the present collapsed fraction. This 
quantity is of interest, for example, since, by combining it with the observed luminosity 
density of the universe, we can get a handle on the average mass-to-light ratio of the universe, 
and the amount of dark matter. The complementary quantity, the uncollapsed fraction, is 
of interest, too, since it determines the amount of matter left behind as the intergalactic 
medium, observable in absorption and by its possible contributions to background radiation. 
A knowledge of the amount of matter left uncollapsed is also necessary in order to interpret 
observations of gravitational lensing of distant sources by large-scale structure. In addition, 
as we shall see, the dependence of the asymptotic collapsed fraction on the equation of state 

* This is a subjective notion, since there is no precise definition of the freeze-out epoch. For a flat, universe with positive cosmo- 
logical constant, for example, spherical density fluctuations must have fractional overdensity S = (p — p) / p > (729pv /50Qp)^/^ 
in order to undergo gravitational collapse, where pv is the vacuum energy density and p is the mean matter density (Weinberg 
1987). In this case, once p drops to a value of the order of pv or less, only density enhancements which are already nonlinear 
will remain gravitationally bound. As such, the "freeze-out" epoch corresponds roughly to the time when p ^ py. Recent 
estimates from measurements of distant Type la SNe, however, suggest values which, if interpreted in terms of this model, are 
closer to p Pv/2 (Garnavich et al. 1998a; Perlmutter et al. 1998), so "freeze-out" began in the past for this model. 
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The asymptotic collapsed fraction in an eternal universe 7 

of the background universe will imply that theoretical tools, such as the Press-Schechter 
approximation, require adjustment in order to take proper account of the effect of "freeze- 
out" on the rate of cosmic structure formation. 

In this paper, we compute the asymptotic collapsed fraction for unbound universes, 
using an analytical model involving spherical top-hat density perturbations surrounded by 
shells of compensating underdensity, applied statistically to the case of Gaussian random 
noise density fluctuations, a model introduced by MSW for the particular case of a flat 
universe with a cosmological constant. We consider a generic cosmological model with 2 
components, a nonrelativistic component whose mean energy density varies as p oc a~^, 
where a is the FRW scale factor, and a uniform, nonclumping component whose energy 
density varies as px a~", where n is non-negative. In terms of the equations of state for 
these two components, we can write this as Pi = WiPiC^, where pi and pi are the mean energy 
density and pressure contributed by component i. For the nonrelativistic matter component, 
w = 0, while for component X, — 1 < w < is the physically allowed range in models in 
which the universe had a big bang in its past and the energy of component X was not more 
important in the past than that of matter, which corresponds to n = 3{l + w) and the range 
< n < 3. The latter condition is necessary in order to be consistent with observations of 
cosmic structure and the CMB anisotropy today. Special cases of this model include models 
with a cosmological constant (n = 0), domain walls {n = 1), infinite strings (n = 2), massive 
neutrinos (ra = 3), and radiation background (n = 4) (although, as explained above, we shall 
exclude values of n > 3 in our treatment here). This generic model, or similar ones, have 
been discussed previously by many authors (e.g. Fry 1985; Charlton & Turner 1987; Silveira 
& Waga 1994; Martel 1995; Dodelson, Gates, & Turner 1996; Turner & White 1997; Martel 
& Shapiro 1998). Recently, such models have been referred to as "quintessence" models (e.g. 
Caldwell, Dave, & Steinhardt 1998) or as models involving "dark energy."[j] The Friedmann 
equation for this model is 



(i-i]o-f^xo)f-)~'+^^of-r'+a 

Van/ Van/ 



'xo 

^ag/ Vag/ Vag 



(1) 



t We note that in some models, the X component is not entirely nonclumping: For massive neutrinos, for example, the 
assumption that the X-component is nonclumping is a very good approximation only for fluctuations of wavelength smaller 
than the "free-streaming," or "damping," length of the neutrinos and for epochs such that longer wavelength fluctuations are 
still in the linear amplitude phase. 
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where H is the Hubble constant, Q = p/pc, = Px/Pc Pc = SH^/SnG, and subscripts zero 
indicate present values of time-varying quantities. 

In §2, we derive the conditions that the cosmological parameters must satisfy in order 
for the background universe to qualify as an eternal, unbound universe. In §3, we compute 
the critical density contrast 6c, defined as the minimum density contrast a spherical pertur- 
bation must have in order to be bound. In §4, we derive the asymptotic collapse fraction 
/c,oo in an unbound universe, using the model introduced by MSW involving compensated 
spherical top-hat density fluctuations. In §5, we compute /c,oo using the Press-Schechter 
approximation, instead. In §6, we compare the predictions of the two models. As we shall 
see, this comparison points up a fundamental limitation to the validity of the ad hoc, over- 
all correction factor of 2 by which the Press-Schechter integral over positive initial density 
fluctuations is traditionally multiplied so as to recover a total collapsed factor which takes 
account of the accretion of mass initially in underdense regions. In particular, we shall de- 
rive this factor of 2 for the Einstein-de Sitter case, but show that the same factor of 2 in 
the Press-Schechter formula overestimates the asymptotic collapsed fraction for an unbound 
universe. To illustrate the importance of these results for currently viable models of cosmic 
structure formation, we apply our model in §6 to two examples of the Cold Dark Matter 
(CDM) model, with fio = 0.3 and Ho = 70kms-^Mpc~\ the open, matter-dominated model 
and the fiat model with cosmological constant. 

2 CRITERIA FOR AN UNBOUND UNIVERSE 

The Friedmann equation (1) describes the time-evolution of the scale factor a{t). The so- 
lutions of this equation can be grouped into four categories, according to their asymptotic 
behavior at late times. If the derivative d, which is initially positive, remains positive at all 
times, never dropping to zero, then the universe is unbound.|^] This is the case, for instance, 
in a matter-dominated universe with Qq < 1. If, instead, d drops to zero as a ^ oo, then the 
universe is marginally bound. This is the case for the Einstein-de Sitter universe {Qq = 1, 
^xo = 0). If d drops to zero at a finite value a = at, then two situations can occur: If the 
second derivative d is negative at a = at, the universe will turn back and recoUapse. This 
is the case for a matter-dominated universe with fio > 1 • However if both d and d are zero 
at a = at, then the universe asymptotically approaches an equilibrium state with a = at 

t The asymptotic value of a in the Umit a —* oo can be either finite or infinite 
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The asymptotic collapsed fraction in an eternal universe 9 

late times. This is the case of the de Sitter universe with a positive cosmological constant, 
which initially expands, and asymptotically becomes an Einstein static universe. 

To determine in which category a particular model falls, we need to study the properties 
of the Friedmann equation (1). For convenience, we rewrite this equation as 

g{y) = H^yf = {1-^0- ^^xo)l/ + ^^o + ^^xol/'"" , (2) 

where y = a/ao = l/{l + z). Only non- negative values of g are physically allowed. Since y > 
after the big bang, the condition g = is equivalent to ?/ = (or a = 0). The first term in the 
right-hand-side of equation (2) can be either positive or negative, while the last two terms 
cannot be negative.^ If n = 2 or fixo = 0, then the quantity fixo cancels out in equation (2). 
This merely illustrates the fact that a universe with a uniform component whose density 
varies as a(t)~^ (cf. a universe with infinite strings) behaves exactly like a matter-dominated 
universe. Such a universe is bound, marginally bound, or unbound if fio > = or 

^7o < I5 respectively. The case in which Qxo 7^ and n = 3 is exactly the same as that with 
^xo = 0, except that Qq is everywhere replaced by Qq + ^xo- In that case, the universe is 
bound, marginally bound, or unbound according to whether Qq + Qxo > 1? ^0 + ^xo = 0, 
or Qq + ^xo < I5 respectively. Let us now consider cases with Qxo 7^ for which n ^ 2 and 

For 1 — f2o — f^xo ^ 0, the universe cannot be bound. Clearly, if 1 — f^o ~ ^xo > 0, then 
g{y) > for all y, and the universe is unbound for any value of n. If 1 — fig^^xo = 0, then the 
last term in equation (2) will eventually dominate (since we assume n < 3). Two situations 
can then occur. If n > 2, then g{y) grows more slowly than y, implying that y"^ = g{y)/H^y 
decreases as y increases, reaching zero as y — > cxd. This is the case of a marginally bound 
universe. If n < 2, y"^ will eventually increase with y. The universe is then unbound. 

Let us now focus on the case 1 — fio — f2xo < 0. If n > 2, then at small y, g{y) > 0, 
but as y increases, the first term in equation (2) will eventually dominate the other terms, 
giving g{y) < 0. There will therefore be a change of sign of g{y) at some finite value y = yt 
where g{yt) = 0. That corresponds to a bound universe. This leaves the interesting case of 
a universe with 1 — fio — ^^xo < and n < 2. Since the slope of g{y) at early times for 
any n < 2 is (1 — — ^^xo) < 0, while at late times it is (3 — n)Qxoy^~"' > 0, g{y) has 
a minimum at some intermediate value of y. When g{y) is zero at that intermediate value, 

§ We are ignoring the possibiUty that Cxo < 0, as would be the case, for instance, in a universe with a negative cosmological 
constant 
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Figure 1. Schematic plot of function g{y) versus y for three different universes if 1 — Qq — f^xo < a^nd n < 2: an unbound 
universe, a marginally bound universe, and a bound universe. The marginally bound case is characterized by the existence of 
a point yt where g = dg/dy = 0. 

this corresponds to the case of a marginally bound universe. The various possibilities for the 
cases with 1 — — ^xo < and n < 2 are shown in Figure 1. The top curve shows a case 
for which g{y) > for all y, that is, an unbound universe. The bottom curve shows a case 
for which g{y) drops to zero at a finite value of y. In this case, the universe turns back and 
recoUapses. It is therefore bound.|5 The transition between these two cases, a marginally 
bound universe, is illustrated by the middle curve in Figure 1, which is tangent to the y- 
axis. At y = yt, both the function g{y) and its first derivative dg/dy vanish. The condition 
for having a marginally bound universe is, therefore, given by the following simultaneous 
equations, 

(1 - ^]o - ^xo)yt + ^0 + nxoyt"" = o , (3) 

{l-Qo- Qxo) + (3 - n)Q^oyt"" = . (4) 

We can solve equation (4) for yt, and substitute this yt into equation (3). We get, after some 
algebra, 

^ At large y, the function g{y) becomes positive again, indicating that there are possible solutions for y large. These are 
"catenary universes," sometimes referred as "no big bang solutions." In such models, the universe contracts from an infinite 
radius, turns back, and reexpands forever. These solutions are not considered to be physically interesting. 
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The asymptotic collapsed fraction in an eternal universe 11 

6 — n J \2 — nj 

We can easily check some limiting cases. For a matter-dominated universe (fixo = 0), equa- 
tion (5) gives r2o = 1 as the condition for a marginally bound universe, as expected. For a 

universe with a nonzero cosmological constant (n = 0), equation (5) reduces to 
27 

(Q, + \, + i)^ = -\,nl (6) 

where we have replaced Qxo by Aq. This is actually a well-known result (see, for instance, 
Glanfield 1966; Felten & Isaacman 1986; Martel 1990). 



3 THE CRITICAL DENSITY CONTRAST 

Consider, at some initial redshift ^ 1, a spherical perturbation of density pi=pi{l + Si) in 
an otherwise uniform background of density pi. Let us focus on positive density perturba- 
tions {6i > 0). Clearly, if the background universe is bound or marginally bound, then the 
perturbation is bound. However, if the background universe is unbound, then the perturba- 
tion can be either bound or unbound depending upon the value of the initial density contrast 
Our goal in this section is to derive the critical density contrast 6i^c, which is defined as 
the minimum value of 6i for which the perturbation is bound. To compute 6i^c, we make use 
of the Birkhoff theorem, which implies that a uniform, spherically symmetric perturbation 
in an otherwise smooth Friedmann universe evolves like a separate Friedmann universe with 
the same mean energy density and equation of state as the perturbation.[tt] Pursuing this 
analogy, a perturbation with 6i > 6i^c behaves like a bound universe, a perturbation with 
6i < 6i^c behaves like an unbound universe, and a perturbation with 6i = 6i^c behaves like a 
marginally bound universe. We can then use the results of the previous section to compute 

First, we need to derive expressions for the "effective cosmological parameters" of the 
perturbation. Notice first that an overdense perturbation has been decelerating relative to 
the background between the big bang and the initial redshift Zi. Hence, at 2; = Zi, the 
perturbation is expanding with an "effective Hubble constant" H- which is smaller than 
the Hubble constant Hi of the background universe. Assuming that the redshift Zi is small 
enough for linear theory to be accurate and for the universe to resemble an Einstein-de Sitter 

II Note: For a nonuniform spherically symmetric perturbation, every spherical mass shell evolves as it would in a universe with 
the same mean energy density and equation of state as that of the average of the sphere bounded by that shell. 
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universe (fij ^ 1, Qxi 1), but late enough to allow us to neglect the linear decaying mode, 
we can easily compute the relationship between H'^ and 6i, 



(7) 

(see, for instance, Lahav et al. 1991). The effective density parameters of the perturbation 
are then given by 

SnGpi 87rGpi{l + 5i) ^(1 + 5.,) 



' 3iff 3/72(1 _ 5./3)2 (l-5,/3)2' ^''> 

" 3iJf " 3/f2(i _ 5^/3)2 - (1 _ s^/s)2 ■ 
Next, we need to find combinations of Q'^ and Q'xi ^^at correspond to "effective" marginally 
bound universes. For the cases for which n < 2, this condition is given by equation (5). We 
now replace fio and fixo by fi^ and fi^j in equation (5)Q and replace Si by 5i^c- This equation 
becomes 



3 — n 



- 3 — 71 


ra(l + (5,Ji 




2-n 



n 



Xi 



(10) 



Since 6i^c ^ 1, we can expand this expression in powers of 6i^c and keep only leading terms. 
We can then simplify this expression further by using the approximation ^ 1. Equa- 
tion (10) reduces to 
fi, + fixi-l + 55i,e/3l3-" nxi 



n 



(2 -n)2-" 



(11) 



6,. 



Xi 



(12) 



Notice that we had to keep the term Qi in the left hand side because of the presence of the 
term —1, and that we cannot expand the left hand side in powers of 6i^c because the quantity 
Qi + Qxi — 1 might be as small as 6i^c- We now solve this equation for 6i^c, and get 
(3 - n)r]]/f 

(2 - ^)(2-n)/{3-n) 

This gives the critical density contrast as a function of the initial density parameters Qi 
and Qxi- We can reexpress it as a function of the present density parameters Qo and 
Qxo and the initial redshift, as follows: The initial density parameters are given by Qi = 
8ivGp,/3Hf = 8rrGpo{l+z,y{Ho/H,f/3H^ = Q^{l+Zif{H^/ H,f and Qx^ = 8irGpx/3Hf = 
8'KGpxo{l + ZiT{Ho/Hif/3Hl = Qxo{l + ZiY{H^/ H if . The ratio {Ho/Hif is given directly 
by equation (1) (with ao/cti = 1 + Zi). We substitute these expressions into equation (12), 
and, using the fact that Zi ^ 1, we keep only the leading terms in (1 + Zi)~^ . Equation (12) 
reduces to 



That equation was derived using the present values of the density parameters, but it is of course valid at any epoch. 
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3 [ (3-n) ^ 1^x0 



) 



/(3-n) 



(13) 



5(1 + 2,) [(2-n)(2-n)/{3-n) ^ 



For the particular cases of a matter-dominated universe {flxo = 0) or a flat universe with a 
nonzero cosmological constant (n = 0, 1 — flo — fixo = 0), we recover the results derived by 
Weinberg (1987) and Martel (1994, eqs. [7] and [8]). 

For cases in the range 2 < n < 3, the condition for a marginally bound universe is 
1 — Qo — ^xo = 0. We substitute equations (8) and (9) into this expression, make the same 
approximations as above, and get 



In terms of the present density parameters and the initial redshift, this expression reduces 

to 



The case n = 3 differs from all others in that the energy density of the X component 
does not diminish relative to that of the ordinary matter component as we go back in time. 
As such, we are never free to assume that the early behavior of the top-hat is the same as 
it would be in the absence of the X component. We shall, therefore, for simphcity, exclude 
this case n — 3 from further consideration here. 

4 THE ASYMPTOTIC COLLAPSED FRACTION 

Our goal is to compute the fraction of the matter in the universe that will eventually end 
up inside collapsed objects (the asymptotic collapsed fraction). Clearly, this question only 
makes sense in unbound or marginally bound universes. In general, the answer depends 
upon the mass scale of the collapsed objects being considered. For cosmological models with 
Gaussian random noise initial conditions (the usual assumption), the density contrast S{X) 
for fluctuations of comoving length scale A is of order [k^P{k)Y^'^, where k — 27r/A is the 
wavenumber, and P{k) is the power spectrum. For a model such as Cold Dark Matter, for 
instance, the power spectrum decreases more slowly than k"^ at large k. Thus the density 
contrast diverges at small scale. Normally, we eliminate small-scale perturbations from the 
calculation by filtering the power spectrum at the mass scale of interest, typically the mass 
required to form a galaxy. The density fiuctuations at that scale have a variance cr^ given by 




(14) 




(15) 




(16) 
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where ll^ is a window function, and R is the comoving radius of a sphere enclosing a mass in 
the unperturbed density field which is equal to the mass scale of interest. Assuming that the 
initial conditions are Gaussian, the fluctuation distribution for positive values of 5 is given 
by 

m) = 47^e-^V2.^ . (17) 

Our problem consists of computing the asymptotic collapsed fraction involving initially 
positive density fluctuations of mass equal to that contained on average by a sphere of co- 
moving radius R, together with the additional mass which eventually accretes onto these 
positive density fluctuations from initially underdense regions, starting from initial condi- 
tions described by equations (16) and (17). In this section, we consider the analytical model 
introduced by MSW. In the next section, we will consider the well-known Press- Schechter 
approximation, instead. 

Consider, at some early time ti, a spherical, top-hat matter-density fluctuation of volume 
V and density contrast 6i, surrounded by a compensating shell of volume U and negative 
density contrast, such that the average density contrast of the system top-hat + shell van- 
ishes. This model is parametrized by the shape parameter s = V/U. If Si > Si^c, the top-hat 
core will collapse. Furthermore, a fraction of the matter located outside the top-hat, inside 
the shell, initially occupying a volume U' < U, will be accreted by the top-hat. Since the 
density is nearly uniform at early times, the asymptotic collapsed mass fraction of this sys- 
tem is simply {V + U') / {V + U). We now approximate the initial conditions for the whole 
universe as an ensemble of these compensated top-hat perturbations, with a distribution 
of top-hat core positive density fluctuations given by equation (17), and we neglect the in- 
teraction between perturbations. As discussed in MSW, the value of s = corresponds to 
the limit in which each positive fluctuation is isolated, surrounded by an inflnite volume 
of compensating underdensity (at a total density inflnitessimally below the mean value p). 
For a flat universe with nonzero cosmological constant, this case was treated by Weinberg 
(1996). The case s = oo corresponds to the limit of "no infall" in which the additional mass 
associated with the compensating underdense volume V is negligible compared with that of 
the initial top-hat. This case was considered for the flat universe with Aq 7^ by Weinberg 
(1987). If s = 1, however, the volume occupied by every positive fluctuation is surrounded 
by an equal volume of compensating negative density fluctuation. This is the case most 
relevant to the problem at hand, involving a Gaussian-random distribution of linear density 
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fluctuations, since the latter ensures that the volumes initially occupied by positive and 
negative density fluctuations of equal amplitude are exactly equal. The full range of values 
of s, < s < oo, was treated by MSW for the flat universe with Aq 7^ 0, with a special 
focus on s = 1 as the case corresponding to Gaussian-random noise initial conditions. The 
insensitivity of the results for the anthropic probability calculations presented there to the 
value assumed for s suggests that the relative amount of total collapsed fraction in universes 
with different values of pv may not be sensitive to the crudeness of the treatment of the 
effect of one fluctuation on another. However, we will also present results here for the full 
range of values of s, while noting that the value s = 1 is the most relevant to the case at 
hand of Gaussian random density fluctuations. 

Under these assumptions, the asymptotic collapsed fraction for the whole universe is 
given by 

where cr, is the value of a at time ti (MSW). For bound and marginally bound universes 
(including, in particular, the Einstein-de Sitter universe), 5i^c = 0, and equation (18) reduces 
trivially to /c,oo = 1 ior all values of s. Hence, the MSW model predicts that, in an Einstein- 
de Sitter universe, all the matter will eventually end up in collapsed objects. For unbound 
universes, we change variables from 5 to x = S'^/2af. Equation (18) reduces to 

f r '"^^"^ (19) 

where 

This equation shows that the collapsed fraction /c^co is unity only when P — 0, which requires 
5i^c — 0. However, in an unbound universe, Si^c is always positive. Hence, according to the 
MSW model, the collapsed fraction in an unbound universe is always less than unity. Notice 
that the dependence upon the cosmological parameters is entirely contained in the parameter 
p. For any cosmological model, we can compute (Tj using equation (16) and Si^c using either 
equation (13) or (15). Since ai oc {1 + Zi)~^ at large Zi for any universe with n < 3, the 
dependence on Zi cancels out in the calculation of /3, as it should: The asymptotic collapsed 
fraction should not depend upon the initial epoch chosen for the calculation. 

The size of the asymptotic collapse parameter /3 determines not only how large or small 
the collapsed fraction is but also how important the increase of collapsed fraction is due to 
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accretion from the surrounding underdense regions. For small values of /9, the asymptotic 
collapsed fraction is close to unity because both the typical positive initial density fluctuation 
and its fair share of the matter in surrounding regions of compensating underdensity collapse 
out before the effects of "freeze-out" suppress fluctuation growth. Hence, in this hmit of 
small /3, "freeze-out" is unimportant and the results resemble that for an Einstein-de Sitter 
universe. For values oi P ^ 1, however, the typical collapse occurs after "freeze-out" has 
begun to limit the growth of density fluctuations. The large f3 limit, in fact, is that in 
which only a rare, much-highcr-than-average, positive density fluctuation is able to collapse 
out of the background before "freeze-out" prevents it, and very little of the compensating 
underdense matter condenses out along with it. For this large /3 limit, equation (19) can be 
shown to reduce to the following simple formula (see Appendix A), 



5 THE ASYMPTOTIC LIMIT OF THE PRESS-SCHECHTER 
APPROXIMATION 

In the Press- Schechter approximation (Press & Schechter 1974; henceforth, "PS"), the col- 
lapsed fraction at time t is estimated as follows: Consider a spherical top-hat perturbation 
with an initial linear density contrast 5i chosen such that this perturbation collapses pre- 
cisely at time t. The density contrast of that perturbation is infinite at time t. However, if 
we estimate the density contrast at that epoch using linear perturbation theory, we obtain 
instead a finite value 6 = Ac, because linear theory underestimates the growth of positive 
fluctuations. The value of Ac is usually taken to be (3/5)(37r/2)^/^ = 1.6865, though this 
result is strictly correct only for the Einstein-de Sitter universe (cf. Shapiro, Martel, & Ihev 
1999, and references therein). A larger perturbation would collapse earlier, and linear the- 
ory would predict that its density contrast exceeds Ac at time t. To compute the collapsed 
fraction at time t, we simply need to integrate over all perturbations whose density con- 
trast predicted by linear theory would exceed Ac at time t, using the distribution given by 
equation (17). The resulting expression, after multiplication by a factor of "2" to correct for 
the fact that half the mass was initially in underdense regions outside the positive density 
fluctuations, is 




(21) 




(22) 
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The introduction of this ad hoc correction factor of "2" in equation (22) is based on 
some assumption about the amount of matter located in unbound regions, either under- 
dense or overdense, which is destined to be accreted onto collapsed perturbations. Consider, 
for instance, the case of an Einstein-de Sitter universe at late time. The critical density 
contrast distinguishing a bound from an unbound density fluctuation is zero, and therefore 
all overdense perturbations are bound and will eventually collapse. Since, for Gaussian per- 
turbations, the overdense regions initially contain only half the mass of the universe, the 
asymptotic collapsed fraction, without taking accretion into account, would be /c,oo = 1/2. 
However, it can easily be shown that in an Einstein-de Sitter universe, all matter in the uni- 
verse will eventually end up inside bound objects. Hence, for this particular case, the proper 
way to handle accretion is to multiply the collapsed fraction by a factor of 2. Equation (22) 
is derived by assuming that this factor of 2 is valid, not only for the asymptotic limit of the 
Einstein-de Sitter universe, but for all universes and at all epochs. Hence, the PS approx- 
imation assumes that the total mass accreted by collapsed positive density fluctuations is 
instantaneously equal to the total mass of these collapsed objects themselves. 

What is the asymptotic collapsed fraction according to this PS approximation? We now 
change variables from 6 to x — P jla^ . Equation (22) reduces to 

fPS -J_r ^ f23) 
where 

To compute the asymptotic collapsed fraction, /^^, we need to take the limit of equations 
(23) and (24) as i ^ oo. Consider a bound spherical perturbation, with the values of its 
initial density contrast 6i at initial time tj chosen so that it collapses at t = cxo. Call this 
value of Si, 5i^oo- By definition, the quantity Ac at t = oo is given by 

Ac(oo) = 5,00 ^4?T ' (25) 
where 5+(t) is the linear growing mode. Since this spherical perturbation collapses at t = 
oo, the initial density contrast Si^oo must be equal to the critical density contrast Si^c- If 
Si,oo was less than Si^c the perturbation would not collapse at all, while if it was greater, 
the perturbation would collapse at a finite time. We can therefore replace 5j_oo by Si^c in 
equation (25). Finally, we notice that the quantity a also evolves according to linear theory. 
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a(oo) = a, +\ / . (26 
Combining these results, we get 

(see eq. [20]). Hence, the PS /3ps parameter reduces to the MSW {3 parameter in the hmit 
t — > oo. Now, comparing equations (19) and (23), we see immediately that these equations are 
identical in the limit s — > oo. Notice that the product 5i,c(^+(oo) takes the undetermined form 
O-oo in the case of an Einstein-de Sitter universe. In this case, the quantity 5i_c<^+(oo) / {U) is 
equal to (3/5)(37r/2)^/^, or 1.6865, at all times. In the Einstein-de Sitter case, (3 = /3ps(cx3) = 
0, and equations (19) and (23) arc the same for all values of s; the asymptotic collapsed 
fractions in that case arc all equal to unity. 

In the limit of large in equation (23), with /3ps replaced by /9, according to 

equation (27), can be shown to reduce to the following simple formula (see Appendix A): 

fcUP»^)-J^.- (28) 
A comparison of equations (21) and (28) reveals that the asymptotic collapsed fraction 
according to the PS approximation is just a factor of (s + l)/s times /c,oo according to the 
MSW model, in the limit of large (3. 



6 DISCUSSION AND CONCLUSION 

Our analjd;ical result in equations (19) and (20) for the asymptotic collapsed fraction in an 
eternal universe can be evaluated for any background universe which satisfies the conditions 
given in §2 which identify it as an unbound universe. Wc need only specify the background 
universe and the power spectrum of primordial density fluctuations, in order to evaluate 
13. Before we do this for a few illustrative cases, however, it is instructive to evaluate the 
asymptotic collapsed fraction /c,oo in general as a function of (3 and s, and compare /c,oo to 
the prediction of the PS approximation, fctoi according to equations (23) and (27). 

We have shown above that the asymptotic collapsed fraction predicted for an eternal 
universe by the PS approximation differs from that predicted here by the spherical model of 
MSW (as generalized to other background universe cases) for s = 1, the value of the shape 
parameter appropriate for Gaussian random initial density fluctuations, with the exception 
of the Einstein-de Sitter universe, for which = /c,oo = 1- For any eternal universe other 
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than Einstein-de Sitter, in fact, the two approaches predict the same asymptotic collapsed 
fraction only if s = oo, instead. The fact that the two approaches generally predict different 
asymptotic collapsed fractions for s = 1 is not surprising, since the PS approximation never 
concerns itself with the fraction of matter which is inside some gravitationally bound region 
and is, hence, fated to collapse out, as the MSW model explicitly does. Instead, the PS 
approximation assumes that, as long as the matter is located within a region of average 
density which is high enough to make it collapse according to the spherical top-hat model, it 
will not only collapse but will also take with it an equal share of the matter outside this region 
which was not initially overdense. This latter assumption is not correct if the under dense 
matter is not all gravitationally bound to some overdense matter. What is perhaps more 
surprising than this disagreement between the two approaches for s = 1 is the fact that they 
do agree for all models if s = oo. 

The fact that in the limit s — oo the MSW model reduces to the asymptotic limit of the 
PS approximation is significant, because the two models are based on different assumptions. 
In the case of the PS approximation, a factor of 2 is introduced to take accretion into 
account. In the MSW model, the limit s — > oo corresponds to perturbations surrounded 
by underdense shells of negligible volume and mass. In this limit, there is essentially no 
accretion. However, the volume filling factor of overdense regions, which is 1/2 in the PS 
approximation, approaches unity in the limit s — > oo for the MSW model, resulting once 
again in a factor of 2 in the expression for the collapsed fraction, but for a different reason. 

We have computed the collapsed fraction predicted by the MSW model as a function 
of the parameter (3, for various values of s, by numerically evaluating equation (19). The 
results are plotted in Figure 2. In addition, the analytical expression in equation (21) which 
is vahd in the large P limit is plotted in Figure 2 for the case s — 1. The analytical expression 
provides an excellent fit to the exact results for the important case of s = 1, not only for 
large (5, but for all j3 ^ 1. The error even at /3 = 1, for example, is only 15%, while at 
P — 5, the error is reduced to 4.5%. For comparison, we also show the prediction of the PS 
approximation, according to equations (23) and (27). This curve is identical to the curve 
for /c,oo for the case s = oo. The point P — 0, fc,oo — 1 corresponds to the Einstein- 
de Sitter universe. As we see, all curves go through this point, indicating that the MSW 
model predicts the correct asymptotic limit in this case, for any value of s. The s = 1 case is 
particularly important, since it is the only one which gives equal filling factors to overdense 
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Figure 2. Asymptotic collapsed fraction /c,oo versus ^8, calculated using the MSW model for various values of the shape 
parameter s (solid curves). The curve for s = oo is identical to ff^, the asymptotic limit (i oo) of the Press-Schechter 
approximation. Also plotted is the simple algebraic formula of the MSW model in equation (21), derived for the large /3 limit, 
for the case s = 1 (dashed curve). 

and underdense perturbations, a requirement for describing realistic Gaussian random initial 
conditions. 

Figure 2 shows that for finite values of s and for /3 > 0, the asymptotic collapsed fraction 
/c,oo predicted by the MSW model is always less than the asymptotic collapsed fraction 
ffoo predicted by the PS approximation. For s < 1, this is not surprising, since in this 
limit the filling factor of the overdense regions is below the value of 1/2 assumed by the 
PS approximation. However, if s > 1, then the filling factor of overdense regions exceeds 
1/2, indicating that the bound perturbations contain more mass in the MSW model than 
in the PS approximation. In spite of this, we still have /c,oo < fc,^- This is caused by 
their different treatments of accretion. The MSW model includes a detailed calculation of 
the amount of matter accreted by a spherical top-hat, while the PS approximation simply 
assumes that the accreted mass equals the initially overdense mass, for all cosmological 
models. Figure 2 suggests that this approximation can be quite crude in some situations 
and greatly overestimate the amount of matter actually accreted. 

To estimate this effect, we have computed, for the MSW model, the "accretion factor," 
Face, defined as the ratio of the total asymptotic collapsed fraction /c,oo divided by the 
asymptotic collapsed fraction /*^ that we would obtain if accretion were neglected. (This 
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100 




Figure 3. Accretion factor Face [= /c,oo (accretion included) //c,oo (no accretion)] versus /3 for the MSW model, for various 
values of the parameter s. The dashed line indicates the value of 2 which is used in the Press-Schechter approximation. 

factor Face is 2 for the PS approximation) . We can easily compute /*^ by going back to the 
derivation of MSW and dropping the term in equation (19) which represents the accreted 
matter. The resulting expression is 



i)/; 



Hence, the accretion factor is 



s + 1 



fc,c 



(29) 



(30) 



s / /c,oo(oo,/3) ■ 

In the large /3 limit, equations (21) and (30) indicate that Fa^{s,P ^ 1) — 1; in this hmit, 
none of the matter in the compensating underdense regions is able to condense out. 

In Figure 3, we plot this accretion factor Fg^c as a function of /3, for various values of s. 
The factor = 2 for the PS approximation is indicated by the dashed line. For the MSW 
model, the accretion factor depends mostly on the amount of matter available in the shell 
surrounding the top-hat, which goes to zero in the limit s — > oo and to infinity in the limit 
s — > 0. For the interesting case s = 1 (underdense and overdense regions with equal filling 
factors), we recover the PS limit Face = 2 at small /3, but the value departs rapidly from 2 
at larger f3. At (3 = 1, for example, the accretion factor drops to 1.125, indicating that the 
PS approximation overestimates the amount of matter being accreted by a factor of 8! 

To demonstrate the importance of this effect for actual cosmological models, we consider 

© 0000 RAS, MNRAS 000, 000-000 



22 

two variations of the Cold Dark Matter (CDM) model: (a) open, matter-dominated CDM 
{Qxo = 0), and (b) flat CDM with nonzero cosmological constant (fixo = Aq = 1— f^o, = 0), 
both with an untilted primordial Harrison-Zel'dovich power spectrum]^ The primordial 
density fluctuation power spectrum for this model, consistent with the standard inflationary 
cosmology and the measured anisotropy of the cosmic microwave background according 
to the COBE DMR experiment, is described in great detail in Bunn & White (1997, and 
references therein). In the absence of tilt, this power spectrum (extrapolated to the present 
according to linear theory) is given by 



P{k) = 27:'[^J 6lk-Tl^M- (31) 
where c is the speed of light and Tcdm is the transfer function, given by 
TcDM(g) = ^-^^^^Y^p^l^ + 3-89g + (le.lg)^ + (5.46g)3 + (6.71g)r'/' (32) 
(Bardeen et al. 1986), with q defined by 

q = (-^)a-n^ohr'Ql,, (33) 



VMpc~\ 

a = ar''''/"''a2-(^''/^°)' , (34) 

ai = (46.9fio/i')°-'™[l + (32.ir]o/i')"°-'''] , (35) 

a2 = (12.0fio/^')°-'''[l + (45.0fio/^')~°-'''] (36) 

(Hu & Sugiyama 1996, eqs. [D-28] and [E-12]), where ^2;, is the density parameter of the 
baryons, and 02. 7 is the temperature of the cosmic microwave background in units of 2.7K. 
The quantity 6h is given by 

Sh = \ 37 

\l.94xlO-5fio°-«°^''^^«, Ao = 1 - 1^0, no tih; ^ ^ 

Once the power spectrum is specified, we can compute the variance cr^ of the present 
density contrast (i.e. as extrapolated to the present using linear theory) as a function of the 
comoving length scale or, equivalently, mass scale over which the density field is smoothed, 
using equation (16). We then compute the parameter /3 using equation (20), where 6i^c is 
given by either equation (13) or (15), and (Xj = (r6+{zi)/6^{0). After some algebra, we get 



q r , , , 1 /Q , ^ .12 

(3= 

where the function ?7(f2o, Aq, z) is defined by 



4^0/ 



(38) 



tt The exponent of the primordial power spectrum, which is unity in the absence of tilt, is usually designated by the letter n. 
It should not be confused by the exponent n used in this paper, which is introduced in equation (1) 
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77(Qo,Ao,^) = (l + ^)^ (39) 

(MSW). In the limit z ^ 1, which wc assume here, the function rj becomes independent of 
z. For flat models (f^o + Aq = 1), MSW derived the following expression: 

5/6 r rXo/Qo dw 1 

(40) 



Jo 



6A 

v{n„ i-no,z » 1) = ^^/^^^ ^ ^^3/, 

For matter-dominated models, we can easily compute the function 77 using the expressions 
given in Peebles (1980). For open models, we get 



(o n ^^1^ 2(1 -Op) , 31^0 , 3Qo , 
r;(i2o,0,^ > 1) = — -< 1 + + — , , In 



1 0^)1/2- -1 



(41) 



5Qo I i-^^o (i-nof/^ 

The fraction of matter eventually collapsed into objects created by positive density fluc- 
tuations of mass greater than or equal to some mass M is entirely specified by the parameter 
f3 evaluated for this mass scale as the density field filter mass. Once (3 is known, we can com- 
pute the asymptotic collapsed mass fractions /c,oo and using equations (19) and (23), 
respectively. We consider models with Hq — 70kms~^Mpc~^ {h — 0.7), = 0.015/i~^ 
(Copi, Schramm, & Turner 1995), and ©2.7 = 1- We have computed cr^ using equation (16) 
with a top-hat window function, 

W(kR)^ (sin kR - kR cos kR) . (42) 

[kR)'^ 

In Figure 4, we plot the variation of the asymptotic collapse parameter /3 with the 
filter mass M [which corresponds to the length scale R in equation (42) according to 
M = 47ri?Vc^o/3, or M/Mq = 1.163 x lO^^^^^^/i^l^o] for two cases of interest: (a) open, 
matter-dominated, flo = 0.3, and (b) flat with cosmological constant, flo — 0.3 — 1 — Aq. 
The value /3 = 1 for these two cases corresponds to the mass scales M/Mq — 3.651 x 10^^ 
(open) and 5.778 x 10^^ (flat), respectively. For Hq — 70kms~^Mpc~^ and density parame- 
ter fio = 0.3 (assuming the shape parameter s = 1, as required for Gaussian random noise 
density fluctuations), the open, matter-dominated CDM model and the flat CDM model 
with nonzero cosmological constant yield mass fractions asymptotically collapsed into ob- 
jects created by positive density fluctuations of mass greater than or equal to the galaxy 
cluster mass-scale IO^^Mq of 0.0361 and 0.0562, respectively. These values of the asymptotic 
collapsed fraction are only 55% of the values determined by the Press- Schechter approxima- 
tion. These results have implications for the use of the latter approximation to compare the 
observed space density of X-ray clusters today with that predicted by cosmological models. 

We have also calculated the asymptotic collapsed fractions /c,oo and as a function 
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log ( M/Mq ) 



Figure 4. Asymptotic collapse parameter /3 versus the filter mass scale M for two CO-B-E-normalized CDM models of interest 
(with f2t,h^ = 0.015, Hq = 70 km s~^Mpc~^ , and a Harrison-Zel'dovich power spectrum): (a) open, matter-dominated (f2o = 
0.3, Aq = 0) and (b) flat, with cosmological constant (f2o = 0.3, Ao = 0.7). 

of Qq (assuming s = 1), for four different filter mass scales M/Mq = 10®, 10^, 10^^, and 
10^^ (notice that the length scale R corresponding to a given mass scale varies with Qo). 
The results are shown in Figure 5. In Figure 5a, /c,oo and /f^ are each plotted separately, 
while in Figure 5b, we plot the ratio /f^//c,oo to demonstrate the extent to which the 
Press-Schechter approximation overestimates the collapsed fraction, especially for cluster 
mass objects and above. 
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APPENDIX A: THE LARGE p LIMIT 
Al The MSW Model 

The asymptotic collapsed fraction according to the MSW model is given by 

s p e-^dx 

If we change variables using x = + w), then equation (Al) reduces to 

" 7rV2 io + 1)1/2 + 1- 
In the limit 1//? ^ 1, we can always find a number a such that 1/ (3 <^ a <^ 1. Since af3 ^ 1, 
we can truncate the integral in equation (A2) ai w = a, because the exponential e~^^ is 
negligible for larger values of w. Hence 

7rV2 io s(w + 1)1/2 + r ^^3) 
Since a <^ 1, the integration variable w is always much smaller than unity, and we can 
replace w + 1 by 1 in the denominator. The resulting integral yields 
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The asymptotic collapsed fraction in an eternal universe 

s \ e-^(l - e"^") 



Since /3a ^ 1, the term e^'^" is negligible. The final expression is 



fcAP » 1) 



.s + lj (7r/5)V2 ■ 
A2 The PS Approximation 

The asymptotic collapsed fraction according to the PS approximation is given by 

A change of variables to w — x^^^ allows us to rewrite equation (A6) as follows: 

2 f°° 

For large /3, 

Combining equations (A7) and (A8), we find 
& » 1) 



in(3) 



1/2 • 
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